Asymptotics of Selberg-like integrals by lattice path counting 
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' We obtain explicit expressions for positive integer moments of the probability density 

, of eigenvalues of the Jacobi and Laguerre random matrix ensembles, in the asymptotic 

^SJ ' regime of large dimension. These densities are closely related to the Selberg and Selberg- 
like multidimensional integrals. Our method of solution is combinatorial: it consists in 

' (— I ' ' the enumeration of certain classes of lattice paths associated to the solution of recurrence 

Q^l relations. 

a , 1 Introduction 



Let be a positive integer and C = [0, 1] be an A^-dimensional hypercube. The integral 

N 



I where A(r) = ni<i<j<v(-^« ~ -^i) Vandermonde determinant, was evaluated by 

Selberg. It has found applications in many different areas of mathematics [H El [3] and 
physics [H 13 |6] , in particular in the theory of random matrices 13 |8] . A recent review of 
applications as well as of the history of the field can be found in [9]. Within the context 
of random matrix theory the constant /3 identifies the orthogonal, unitary and symplectic 
^ ■ universality classes and takes value in the set {1,2,4}. The normalized function 

b : N 

P^iT) = -\A{T)fl[TrHl-Tir-' (2) 

is the joint probability density of the eigenvalues of matrices from the Jacobi /3-Ensemble 
[8]. The average value of any function of the eigenvalues f{T) is given by the multiple 
integral 

{fin = [ /(T)P;'T(T)d^r, (3) 

withd^T = ll.dTi. 

One particular application involves quantum electronic transport in mesoscopic struc- 
tures [5]. If there are A^i incoming channels and outgoing channels, the system may be 
described by a x transmission matrix t whose element tij is the probability amplitude 
of transmission from channel j to channel i. For systems with chaotic classical dynamics 
the random matrix approach to the problem consists in assuming the system's unitary S- 
matrix to be a random element of a circular /3-ensemble. This is equivalent to assuming the 
hermitian matrix T = tt^ to be uniformly distributed in the Jacobi /3-Ensemble [TOj with 
N = min{iVi,iV2}, a = f (|iV2 - iVi j + 1) and 7 = 1. 
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In the above context (with 7=1) the average value of quantities hke 

N 

1=1 

can be used to characterize universal statistical properties of quantum chaotic transport. 
Exact results appeared for small n in |1H \T1[ \T3[ [T3] and for general n in \15\ [16] . In this 
last work the present author used a result of Kaneko [2] and Kadell [17] which gives the 
average value of any Schur function of the eigenvalues, an approach that was later taken 
further in [18^ I19j . For large numbers of channels, a generating function for the average 
value of ^ was presented in pO], and an explicit expression appeared in [13] . 

We are concerned with the asymptotic regime 3> 1 of the average value of ([4]) with 
respect to the probability density ([2]). Special cases of this problem were recently discussed 
in [21], by using the method from [16j and taking ^ 1 in the last step. It was then noticed 
that the results had combinatorial interpretations, but no reason for that was provided. The 
authors of [21] also conjectured that in this regime the average has a factorization property, 
i.e. {f{T)g{T)) becomes asymptotically equal to the product {f{T)){g{T)) if the functions 
f{T) and g{T) depend on disjoint sets of variables. If this is true, then the quantity (T") 
becomes indeed the most interesting one to compute. 

We also consider the asymptotic regime ^ 1 of the average value of @ within the 
Laguerre ensemble of random matrices. It that case the eigenvalues belong to [0, 00) and 
have a joint probability density given by 

1 ^ 

where Y{a, (3, e) is a known normalization constant which is given by a Selberg-like integral. 
This also finds an application in the area of quantum chaotic transport \22\ 123] . The 
eigenvalues of the Wigner-Smith time-delay matrix Q = —ihS^^dS/dE (where E is the 
energy) are called proper delay times, T,. The inverse delay times l/Tj are distributed 
according to the Laguerre ensemble with a = (3N/2 and e = —(3N/[2'y) where 7 is the 
classical decay rate of the system. 



2 Statement of results 

We obtain an explicit formula for (T"), in the asymptotic regime, valid for arbitrary values 
of /3 (not restricted to the set {1,2,4}) and for arbitrary parameters a and 7, which are 
allowed to grow linearly with as is sometimes required. Let [x] denote the integer part 
of X and let 

c,= n^ (0) 



n / n + 1 

be the Catalan numbers. We show that, asymptotically, 



1 

n-l / -,\ L 2 



m=0 ^ ^ fc=0 ^ ^ 

where 

A 2a + /3N AAA 

2{a + 'j + pN) 2{a + 'j + pN) 

The above formula generalizes all special cases that have so far been considered. Our 
derivation is essentially combinatorial, and consists of finding a recurrence relation and then 
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turning its solution into the problem of enumerating certain lattice paths. In the course 
of the calculation we prove the factorization conjecture already mentioned, for polynomial 
functions /(T) and g{T). 

For the Laguerre /3-Ensemble ([5]) we show that the factorization conjecture holds as well 
and that the relevant asymptotic average value is given by 

(Tf) = ^2 ^ (''~J\cMA2rA^-'-^"' (9) 

where 

BN 

Ai = ^, A2 = a + Ai, ^3 = ^1+^2. (10) 

We also consider the analogous problem for the distribution of proper delay times. It was 
shown in [24J, by other means, that the solution contains the Schroder numbers; we present 
a combinatorial proof of that. 

We note that the Laguerre and Jacobi ensembles are also widely studied in the context 
of multivariate statistical analysis [H [25] , where they are known as the Wishart distribution 
and the multivariate beta distribution. Moments of the latter were considered in [26\ . 

The paper is organized as follows. In Section 3 we consider the quantum transport case, 
which involves Jacobi Ensembles with 7 = 1. This is the simplest version of the problem. 
The lattice paths involved are Dyck paths, which have only two types of steps. In Section 
4 we turn to the Laguerre Ensembles, which are second in the complexity scale. The lattice 
paths involved are now Motzkin paths, which have three types of steps. In Section 5 we 
consider the distribution of proper delay times in quantum chaotic scattering (closely related 
to Laguerre). The lattice paths that appear in that case are Schroder paths. Finally, in 
Section 6 we come to Jacobi Ensembles with arbitrary 7, which involve lattice paths with 
four types of steps. 



3 Jacobi Ensembles with 7 = 1 

Let /(T) = lY^^^ We will show that {f{T)T^) ^ {f{T)){T^) in the asymptotic regime 
N ^ 1. To that end let us define Q = f {T)T^~^^ Pj^ (T) , with n > Aj. Taking into account 
that 

dP^ (a -I) 1 



dTrr 



p-(T) + l3P^{T)^-^_ 



T —T- 



we can write the derivative of Q with respect to T, 
dQ 



dTrr 



Notice that 



f{T)P^{T) 



SQ 



T. 



n+l 



Ti 



dT„ 



dTrr 



(12) 



(13) 



does not depend on the value of n. We will use this fact to find a recurrence relation. This 
idea goes back to Aomoto [27] . 
The integral of ^ is 



/ ^d^r = (n + a)(/(r)r«) + /3^ 



/ /TO 

\ Tm Ti I 



n+l 

m 



(14) 



But 



//(r)r,r' 



By symmetry. 



+ (iV - m) 



/(r)r, 



-n+l 



T 



m+1 



(15) 



T 



m+l 



Also, in the same vein 



//(T)r, 



n+l 



\ Tlr 



1 / T-in+1 _ 



n 

-5;(/(r)rr'^r^+i). 

a=0 



- 5^ (/(r)T;"rr">. 

a=0 



(16) 
(17) 



(18) 



+ 0(1). 



(19) 



Every term in (jlSp is of the same order in the hmit N oo, and we can thus make the 
approximation 

^ / /(T)r^+^ _ /(r)r^+^ 

This is the crucial step towards the factorization property, because we have just ignored 
the terms of the sum containing the variables that appear in f{T). They are the ones that 
obstruct the factorization for finite N. 

From now on we simply ignore the error terms and keep in mind that the limit N ^ oo 
is implicit. Defining the function 



n-l 



a=l 



we have 



dX, 



(20) 



(21) 



where we have neglected n against a + f3N (remember that a may grow with A^). 

We have seen that the above does not depend on the value of n. In particular, for n = 
it gives {a+ /3N/2){f{T)) and forn = 1 it gives {a+f3N){f{T)Tm). Comparing the formulas 
for general n and n = we obtain 



(/(r)0 = ^2(/(T))-^i(/(W 



m,m+l\ 



where 



^1 



2a + /37V 



(22) 



(23) 



2(a + /3iV)' 2{a + PN)' 

Recurrence relation ()22p proves the factorization conjecture, at least for polynomials. This 
is because the function /(T) does not change when the relation is iterated, so it should be 
clear that it results in (/(r)T^) = {f{T)){T^). 

The use of letters Ai,A2 comes from the quantum transport setting, where if ^ 1 
we obtain Ai ~ Ni/{Ni + A^2)- Also, in that case a is proportional to /3, causing the final 
result to be independent of /3. Notice that if a is held fixed in the limit of large A^ we end 
up with Ai = A2 = 1/2. 



A 



We now turn to our main objective, which is the calculation of the basic quantity (T"). 
We omit the index of the variable since it is irrelevant. Let us define 



n-l 



p„ = ^(r"-'^)(r'^). (24) 



a=l 



We have (T'') = 1, (T^) = A2 and the recurrence relation 

(T") =A2- AiVn, n>2. (25) 

We may write this as 

(T") = A2-^iM„ + 2Ai(r"), n>l, (26) 

where we have defined 

n 

M„ = ^(r"~")(T'^). (27) 

a=0 

It is possible to obtain the ordinary generating function 

F{x) = 1 + ^(r")a;". (28) 



n>l 



Since 



we get the algebraic relation 



[F(x)]2 = l + J^M„x'^ (29) 



n>l 



F = l + ^^ + 2Ai{F-l)-Ai{F^-l), (30) 
1 — X 



which can be solved to give 



F(x) = l + + — (31 

^ ^ 2Ai 2Ai V 1 - X ^ ^ 

The generating function of course provides (T") for any value of n. However, an explicit 
formula for this quantity can be obtained by representing the recurrence relation pictorially 
and turning its solution into a lattice path counting problem. 

Let us start with a brief example. The recurrence relation (I25p gives 

(Tf) = A2- Ai{{TfT2) + {TfTi) + (TiTD). (32) 

The coefficient of Ai contains all ordered partitions of 4 into two parts. In the next step we 
have 



Ao-A^ 



A2i{Ti) + (Tf) + (Tf )) - ^i((riT2r|) + (TiT^n) + (TfTaTg)) . (33) 



Now the coefficient of Af has the ordered partitions of 4 into three parts. The coefficient 
of A1A2 has the ordered partitions -into one part- of the numbers up to three. At every 
step, when the power of Ai is increased the number of parts in the ordered partitions also 
increases; on the other hand, when the power of A2 increases one part of the previous 
partitions is removed. 



[n,4] 



[n,l] 




Figure 1: The recurrence relation ()25p may be mapped into a path-counting problem if we 
associate with Ai and A2 different directions of movement. In the Ai direction we increase 
the number of parts of the ordered partitions, while in the A2 direction we eliminate one 
part. If a path hits the lowest horizontal line it stops. The number of terms is constant 
along falling steps. 



Let [N, m] denote the set of all ordered partitions of into m parts. We can represent 
a general step in the iteration of the recurrence relation as 



[N, m] ^[N,m + 1] 

iA2 iA2 



N-l 



N-1 



[n,m — 1] — 1^ [n, m] 



n=l 



n=l 



Notice the commutativity of the diagram. We may therefore think of Ai and A2 as directions 
in which we can move as we proceed with the calculation. If we further simplify notation 
by writing (A^, m) = [jn=i I'^' "^1 then Figure 1 codifies the recurrence relation. 

We start at [n, 1], which just represents (T"). Iteration of the recurrence relation then 
produces all possible sequences of rising and falling steps that remain above the (-,1) hori- 
zontal line. A path ends if and only if it drops below that line. The final result for (T") will 
therefore be equal to A2 times a polynomial in {—A1A2) of degree n — 1. The coefficient of 
{—AiA2)^ contains two contributions. First, the number of terms in the set [n,p + 1] (be- 
cause moving in the A2 direction does not change the number of terms) . This is the number 
of ordered partitions of n into p + 1 parts and is given by ("p^)- Second, the number of 
different paths leading from [n, 1] to (n — p, 1). Clearly, the relevant paths are Dyck paths: 
lattice paths with steps (1, 1) and (1, —1) that never fall below the x axis. The number of 
such paths containing 2p steps is well known to be the Catalan number, Cp. In conclusion. 



(34) 



This calculation has been sketched previously in jlSj . 



4 Laguerre Ensembles 

Before tackling the general Jacobi ensembles, we first consider the Laguerre case. We now 
have the kernel 



1 ^ 

L^/{T) = y\AfllT,^e-^^', (35) 



1=1 

whose derivative is 



dL n ' C\ I V 1 



^Lf{T) + /3L^'^(T) 5^ ——- - eLf{T). (36) 



The integration domain is now C = [0, oo)^. If we define Q = f {T)T^^ V^'' {T) it is easy to 

see that the value of -^^dTm is always zero, for any value of n. Proceeding analogously 
to what we did to arrive at Eq. ()2ip . we find that 



dQ 



dT„ 



d''T 



+ ^(/(T)z)-'-+i) - 6 (/(r)C+i> = 0. (37) 



The above equation implies that the factorization conjecture holds again, since /(T) is 
not affected by the iteration. Let us define the new variables 

^1 = ^, ^2 = ^ + ^1, vis = ^1+ via. (38) 

Eq. (j37p gives in particular (T) = A2 and (T^) = A^l^). For general n we have 

(r") = A3(T"-i)+vliP„_i, (39) 

where T>n has been defined in Eq. (j24p . Telescoping this equation we obtain the recurrence 
relation 

n-2 

(T") = A^Al-^ +MY, AlVn-k-i- (40) 

fc=0 

In terms of the generating function F{x) = X]n>o(-^")^" elementary to derive the 
simple quadratic algebraic relation 

^^^^M^l±AMF-^^ (41) 
1 - visa; 

The recurrence relation (|40p can also be interpreted in terms of lattice paths. Again 
the constants Ai, A2 and A^ are associated with different directions of movement in the 
plane, which take us from one set of ordered partitions to another. In the A2 direction we 
simply remove one of the parts, as in the previous section. In the Ai direction we increase 
the number of parts by one and decrease the number being partitioned by one. In the A^ 
direction we keep the number of parts constant and decrease the number being partitioned 
by one. If we denote, as in the previous section, by [n, m] the set of all ordered partitions 
of n into m parts, and define the union {N,m) = UnLit""' ""^l' then the iteration of (j40p is 
depicted in Figure 2. 

As a brief example, consider 

(Tf) = vi2^^ + Ai((Tfr2) + (rfr|) + {nri) 

+ ^3(r2r2) + AsiTiTi) + Al{TiT2)). (42) 
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Figure 2: The recurrence relation for the Laguerre ensemble (j40p may be mapped into a 
path-counting problem if we associate with Ai, A2 and ^3 the rising, falling and horizontal 
directions of movement, respectively. Notice that the action of A\ is different than in the 
previous section and in Figure 1. 

We interpret this in the following way: we can take two horizontal steps and one rising step 
to go to the partitions of two into two parts, A\A^l!T\T2)\ we can take one horizontal step 
and one rising step to go to the partitions of three into two parts, ^iA3((TfT2) + (TiTl)); 
we can take only one rising step to reach partitions of four into two parts; finally, we can 
take four horizontal steps and one falling step. In general, from any given point we may go 
up or down, but before we do so we can take any number of horizontal steps. 
With one more iteration we find 

(Tf ) = A2A\^A<,A^{Tl) + A^A^A^ljl) + A<,A\A^{T^) + Al^T^T^T^) 

+ A^A^A^l^Tl) + A2^3^i^3(ri) + ^2AiAi(Ti). (43) 

It is instructive to treat the variables as non-commutative because that makes it easier to 
visualize the paths. So even though A2A'^Ai, 742^3^1^3 and A2AiA^ are all equal, the 
lattice paths to which they correspond are different (notice that the steps should be read 
from right to left). 

The general structure of the calculation of (T" ) is as follows. The last step is always ^2- 
The rest of the path consists in n — 1 steps, which clearly result in Motzkin paths: lattice 
paths with steps (1, 1), (1, —1) and (1, 0) that never fall below the x axis. The number of 
Motzkin paths of length n containing exactly m rising steps is given by 

Mn,m=(^^C^, (44) 

where Cm are the Catalan numbers. In conclusion, 

[^] 

(T") = A2 ^ M^.^^miAiA^TA^-^-^^. (45) 

m=0 

Notice that if limiv^oo f = we have (T") = C„(r)". 
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5 Proper delay times 



It was shown in [22^ [23] that the proper delay times Tj associated with quantum scattering 
by a chaotic cavity are distributed according to 



1 



N 



(46) 



i=l 



with a = —3f3N/2 + 13 — 2 and e = —I3Nt£)/2, where is the cavity's classical dwell time 
and is the number of decay channels. The average value of T" was computed in [Mj by 
integration against the displaced semicircle eigenvalue density, and was found to be related 
to the large Schroder numbers. 

The calculation in this case is quite similar to that in the previous Section. The derivative 
of the kernel is 

a ^ \ - 1 



dTm. 



(47) 



Defining Q = f {T)T^~^^ Z^'"^ (T) , the value of -^^dTm is again zero for any value of n. 
Now we take n > and get 



(a + /3iv)(/(r)0 + ^(/(r)i)-'-+i) + e (/(T)rr^: 



0. 



(48) 



The above equation implies that the factorization conjecture holds again, since f{T) is 
not affected by the iteration. Defining the new variables 



1, A2 = A-s 



2{a + /3N) ^' a + f3N 

Eq. dM]) gives (T) = A2 and, in general, (T") = ^3(r""i) + AiVn- Telescoping gives 



(49) 



n-2 



(50) 



k=0 



The recurrence relation (I50p is only slightly different from (I40p . The difference is that 
steps in the Ai direction no longer change the number being partitioned. The iteration 
of (j50p may be interpreted as in Figure 3: Ai corresponds to a vertical step, while A2 
and A^ correspond to falling and horizontal steps, respectively. The paths involved in the 
calculation of (T") are those going from [n, 1] to [1, 1] without falling below the initial 
horizontal level. Once the step reaches [1, 1] a final A2 step terminates it. The total number 
of A2 steps plus the total number of A^ steps is always equal to n. Since Ai = 1 and 
A2 = A3 = td the value of (T") will be equal to times the number of possible paths. 

The solution to the enumeration problem posed above is nothing but the (large) Schroder 
number In order to see this, suppose we turn every Ai step from vertical to rising (i.e. 
from (0, 1) to (1, 1)) and double every A3 step (i.e. from (1,0) to (2,0)). The resulting path 
will always be a Schroder path of length 2n. The map is bijective, so Rn is the number we 
are looking for. In conclusion, we have 



(51) 



Q 




Figure 3: The recurrence relation (|5U|) associated with proper time delays is mapped into a 
path-counting problem by taking Ai, A2 and A3 to be vertical, falling and horizontal types 
of steps, respectively. These paths are in bijection to Schroder paths. 



6 Jacobi Ensembles 

We now come back to the general Jacobi Ensembles, when the Selberg kernel is 

1 ^ 

P^'^T) = -|A|'3 []7;"-i(l - Ti)^-^ (52) 

i=l 

and the integration domain is again C = [0, 1]^. In this case we have 

dP^ (a -I) ^1 Pr'^{T) 

and the integral of dQ/dTm with Q = f {T)T^^^ Pg''^ (T) is given by 



c dTjn 



a + /3iV ( 1 _ ^ 



+ ^(/(T)D™'™+1) - (7 - 1) (/m^^) . (54) 

It remains independent of re. We have again neglected re compared to a + I3N, and we may 
also approximate 7 — 1 ?a 7 since what matters is limAr_s.oo 7/A^. 
Let us define 

M = — -, A2 = —, r, A-i = Ai+ A2. (55) 

2(a + /3iV + 7)' 2(a + /3iV + 7)' ' ' ^ ^ 

Comparing the values of at n = 1 and re = we have {f{T)Tm) = A2{f{T)). Comparing 
in general re and re — 1 we have 

(/(r)0 = A3(/(r)rr^) + Ai(/(r)[D::!ir^ - d^'"^+']). (56) 

The factorization property clearly continues to hold for general 7. 

Let us now consider (T"). We have (T) = ^2- For general re we may telescope the 
previous equation to obtain 

n-2 

(T") = A2A^-' + (1 - As)Ai J2 A^n-k^i - AiVn, (57) 

k=0 



in 



where T>n has been defined in Eq. ()24p . For the generating function F{x) = ^,j>o(^")^" 
this leads again to a quadratic algebraic relation 

In order to turn the recurrence relation into a lattice path problem it is convenient to 
write it as 

n-2 

(T") = ^2^r' + ^4 5Z ^3^n-fc-l - AiVn, (59) 
k=0 

where A4 = Ai{l — A^) is treated as an independent variable. Comparing this equation 
to the ones obtained in the previous sections, we see that Ai and A2 have the same role 
they had in Section 3, namely A2 removes one part of every partition and Ai increases the 
number of parts by one. ^3 has the same role it had in Section 4: it decreases the number 
being partitioned by one. Finally, A4 now does what Ai did in Section 4, it decreases the 
number being partitioned by one while increasing the number of parts by one. The general 
structure is as depicted in Figure 4. 

We now have lattice paths with four possible steps: (0,1), (1,-1), (1,0) and (1,1) 
corresponding respectively to Ai, A2, A3 and A4. The total displacement in the horizontal 
direction must be n — 1. Let us consider the situations in which exactly n — 1 — m of 
the horizontal moves come from A3 steps. There are ("~ ) possibilities for assigning their 
position. Let us suppose that exactly k of the remaining steps are of type A4. We are 
thus reduced to counting the number of paths with k raising steps, m — k falling steps and 
m — 2k vertical steps. The key to this enumeration is that every one of these paths may be 
obtained from a Dyck path of length 2m — 2k, if we replace m — 2k its raising steps by 
vertical ones. The number of such Dyck paths is Cm-kt ^-nd the number of possibilities to 
single out m — 2k of the raising steps is C"^*^) • 

In conclusion, the final result is given by 



rm+i 

n—1 / 1 \ L 2 



{T-) = A2Y,(''j)Ar'-"^ ^ '')c^^,{-A,A2r'"'{A,A2)'. 



(60) 



Replacing back ^4 by Ai( 1—^43 ) we arrive at the result ([7]). Naturally, when limAr^oo 7/-^ 
we have A3 = 1 and this reduces to 




Figure 4: The recurrence relation for the Jacobi ensemble corresponds to lattice paths where 
Ai, A2, A3 and A4 are respectively associated with vertical, falling, horizontal and rising 
steps. 
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